Entanglement versus energy in quantum spin models 
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We study entanglement properties of all eigenstates of the Heisenberg XXX model, and find that 
the entanglement and mixedness for a pair of nearest-neighbor qubits are completely determined 
by the corresponding eigenenergies. Specifically, the negativity of the eigenenergy implies pairwise 
entanglement. From the relation between entanglement and eigenenergy, we obtain finite-size behav- 
iors of the entanglement. We also study entanglement and mixedness versus energy in the quantum 
Heisenberg XY model. 
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Quantum entanglement lies at the heart of quantum 
mechanics, and can be exploited to accomplish some 
physical tasks such as quantum teleportation [1] . In this 
sense, it can be regarded as a resource, just like energy. 
As pointed out by Osborne and Nielsen [2], the simi- 
larity between entanglement and energy turns out to be 
more than superficial. It is interesting to explore the rela- 
tionship between these two resources, entanglement and 
energy. 

Recently, the study of entanglement properties in 
many-body systems have received much attention [8? - 
19]. Specifically, for the ground state (zero temperature) 
of a ring of N qubits interacting via the antiferromag- 
netic isotropic Heisenberg Hamiltonian, a direct relation 
is established between the concurrence C [20] quantifying 
the two-qubit entanglement and the ground-state energy 
per site eo [10, 21, 22]: 

Co(iV) - max[0, ^£o{N)/N] = max[0, ~eoiN)], (1) 

where C refers to the concurrence for two nearest- 
neighbor qubits, and So is the ground-state energy. For 
a pair of qubits, the entanglement of formation can be 
obtained from the concurrence C 



E, 



1 + Vl - C2 



o/ 



(2) 



where h{x) — —x log2 x — [l — x) log2(l — x). The concur- 
rence itself is a good measure of two-qubit entanglement, 
and we adopt it as our measure of pairwise entanglement. 
For the case of finite temperature, the concurrence for 
two nearest qubits in the thermal state of the ring was 
found to be related to the thermodynamical function, the 
internal energy U via [21, 22] 



C(iV) = max(0,-C//7V), 



(3) 



which connects the microscopic quantity, the concur- 
rence, and the macroscopic quantity, the internal energy. 
Equations (1) and (3) show that the two-qubit entan- 
glement and energies are closely related in the isotropic 
Heisenberg model. 



Having known the relations between entanglement and 
energy for ground states and thermal states of the XXX 
model, further questions arise that what are the entan- 
glement properties of excited states, and what are the re- 
lationships between entanglement and excited-state en- 
ergies. We will address the question in this paper. To 
study entanglement of excited states is not only interest- 
ing itself, but also help to understand the entanglement 
at finite temperatures. We also consider the entangle- 
ment versus energy in the quantum XY model, and show 
that the entanglement exhibits a symmetry, distinct from 
that in the XXX model. 

There exists another concept, the mixedness of a state, 
related to entanglement, is also central in quantum infor- 
mation theory [23]. For instance, Bose and Vedral have 
shown that entangled states become useless for quantum 
teleportation on exceeding a certain degree of mixed- 
ness [24]. For a pure bipartite state, the mixedness of 
one subsystem is equal to that of another one, and can be 
used to quantify bipartite entanglement. In this case, en- 
tanglement and mixedness are equivalent. We will study 
both entanglement and mixedness properties. 

Entanglement and eigenenergy. We consider a physical 
model of a ring of N qubits interacting via the isotropic 
Heisenberg XXX Hamiltonian 



H 



N-l 



i+l 



J S]si , 



(4) 



where Sj,j+i —^ {1 + a^ ■ (Ji+i) is the swap operator be- 
tween qubit i and j , (Tj ~ (aix , <Jiy , Oi^ ) is the vector of 
Pauli matrices, and J is the exchange constant. Posi- 
tive and negative J correspond to the antiferromagnetic 
and ferromagnetic case, respectively. Note that we have 
assumed the periodic boundary condition. 

When an energy level of our system is non-degenerate, 
the corresponding eigenstate is pure. When a /c-th energy 
level is degenerate, we assume that the corresponding 
state is an equal mixture of all eigenstates with energy £k ■ 
Thus, the state correspoding to the fc-th level with degen- 
eracy becomes a mixed other than pure, keeping all sym- 
metries of the Hamiltonian. A degenerate ground state 



is called thermal ground state in the sense that it can 
be obtained from the thermal state exp[— i//(fcBT)]/Z 
by taking the zero- temperature limit [2]. Here, ks is the 
Boltzman constant, T the temperature, and Z the parti- 
tion function. The k-th eigenstate pk can be considered 
as the thermal ground state of the nonlinear Hamiltonian 
H' given by 



H' = {H- £u) 



(5) 



Note that Hamiltonian H' inherits all symmetries of 
Hamiltonian H . 

Due to the translational invariance of the system, all 
nearest- neighbour entanglements are identical. Thus, 
from now on, we consider the entanglement between 
qubit 1 and 2, representing nearest-neighbour pairwise 
entanglement. Another important SU(2) symmetry of 
the Hamiltonian guarantees that the reduced density ma- 
trix p). ' = Tr3^4^... ,N{pk) for fc-th eigenstate is given by 



pr 



f Uk+ \ 

Uk- Zk 

Zk Uk- 

\0 Uk+ J 



the standard basis {|00>, |01), |10>, |11>}, with 



uk± = (1 ± Gr)/4, Zk = Gr/2. 



(6) 



(7) 



Here, G^" 



Tr((TiQ,cr2aPfe), a G {x,y,z} are correlation 
functions which are equal due to the SU(2) symmetry. 

From Hamiltonian H, a simple and useful relation be- 
tween eigenenergy per site e^ and the correlation function 
is obtained as 



ek = £k/N 



1 + 3Gf 



(8) 



Therefore, from Eqs. (6-8), the reduced density matrix 
is completely determined by the eigenenergy e^. As all 
information about pairwise entanglement and mixedness 
of pj.^ are contained in pj.^, the eigenenergy completely 
determines the entanglement and mixedness. Next, we 
give a quantitative relation between the concurrence and 
the eigenenergy. 

From the two-qubit reduced density matrix (6), the 
concurrence for two nearest qubits is obtained as [10] 

Gfc(iV)=2max(0,|zfc|-Wfe+) 

-max(0,|Gr|- Gr/2 -1/2) 

= max(0,-3Gr/2-l/2) 

= max(0,-efc). (9) 

The third equality follows from the inequality jG^ | < 1, 
which is a special case of a more general result that 
I (A) I < 1 for any Hermitian operator A satisfying A^ — I. 
The last equality is obtained by using Eq. (8). Thus, we 
get a simple relation between the pairwise entanglement 



TABLE I: The concurrences of the ground state, first and 
second excited states for the number of qubits from 3 to 16. 



iV 


3 


4 


5 


6 


7 


8 


9 


Co 


0.0 


0.5 


0.2472 


0.4343 


0.3158 


0.4128 


0.3438 


Ci 


0.0 


0.0 


0.0 


0.2060 


0.0161 


0.2821 


0.1525 


C2 


- 


0.0 


0.0 


0.0 


0.0055 


0.1749 


0.1342 


N 


10 


11 


12 


13 


14 


15 


16 


Co 


0.4031 


0.3580 


0.3979 


0.3661 


0.3948 


0.3712 


0.3928 


Ci 


0.3184 


0.2258 


0.3386 


0.2695 


0.3509 


0.2976 


0.3590 



C2 0.2541 0.2100 0.2962 0.2570 0.3212 0.2877 0.3371 



of two nearest qubits and the corresponding eigenenergy. 
A necessary condition for non-zero entanglement is that 
eigenenergy is negative, i.e., the negativity of eigenen- 
ergy implies pairwise entanglement. We also see that the 
ground state exhibits largest entanglement. 

Now we apply the above result to investigate finite-size 
behaviours of the entanglement. Conformal invariance 
theory predicts the finite size behaviour of the ground 
state and the first excited state [25] 



ei{N) - eoiN) K TT^b/N^ 



(10) 



with b being identified with the central charge. From the 
above equation and Eq. (9), we immediately have 



Ci{N)^Ca{N)K-7r^b/N'^ 



(11) 



which is the finite size behaviour of the entanglement. 

In Table I, we give the concurrences of the ground 
states, first and second excited states for the number of 
qubits from 3 to 16. The ground state becomes entangled 
when iV > 4 [21, 22]. However, the first (second) excited 
state becomes entangled when TV > 6 {N > 7). The 
concurrence of the ground state decreases (increases) as 
even N (odd N) increases. However, the concurrence of 
the first excited states increases as even N or odd N in- 
creases. This fact also holds for the second excited states. 
We see that the entanglement properties of the excited 
states are distinct from those of the ground states. 

To clearly display behaviours of the pairwise entangle- 
ment corresponding to different energy levels, in Fig. 1, 
wc plot the concurrences versus energy level index k for 
different A^. It is evident that for a fixed N the con- 
currence decreases monotonically with the increase of k, 
and when k is equal or larger than a threshold value fcth, 
the entanglement vanishes. The threshold values fcth is 
A-dependent, and the larger the number of qubits, the 
larger the threshold value. More and more excited states 
become entangled when the number of qubits increases. 

Next, we study the mixedness properties of all eigen- 
states. The mixedness of a state g can be quantified by 
the linear entropy given by 



E ^1- Ti{g^ 



(12) 
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FIG. 1; The concurrence versus energy level index k for dif- 
ferent number of qubits. 
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FIG. 3: The concurrence versus eigenenergy per site et for 
different number of qubits in the XY model. 



the mixedness of the two-qubit state. For instance, when 
A^ = 6, the ground state, the second excited state, and 
the 10-th excited state are non-degenerate. 

In the isotropic Hcisenberg model described above, 
both entanglement and mixedness of an eigenstate are 
completely determined by the corresponding eigenenergy. 
This result is due to the many symmetries in the model. 
For other models, the relation between entanglement and 
energy may become more complicated. We now consider 
another well-established model, i.e., the Heisenberg XY 
model. The XY Hamiltonian is written as 



200, 400 

k 



H^ 



XY 



J N 



' Ci+ii 



Jiy 



CTj+lj;), (14) 



FIG. 2: The linear entropy versus k for different A'^. 

Then, from Eqs. (6), (7), (8), and (12), the linear entropy 
of the two-qubit state pl^ is obtained as 



Ek{N) 



1 



[4iN)-ek{N) + l] 



(13) 



which is determined solely by the eigenenergy e^ per 
site. Obviously, the linear entropy takes its maximum 
3/4 when efe == 1/2. 

In Fig. 2, we plot the linear entropy versus k for differ- 
ent number of qubits. As k increases, the linear entropy 
monotonically increases and reaches a maximum, then 
monotonically decreases. In contrast to the concurrence, 
the linear entropy takes its maximum in the middle of 
energy levels, other than at the ground-state level. For 
instance, when N = 6, there are 13 levels and the linear 
entropy takes its maximum when k = 6. There are some 
eigenstates which are non-degenerate, and thus pure. In 
this case, the entanglement of between the two qubits and 
the rest quantified by the linear entropy is equivalent to 



where the periodic boundary condition is assumed again. 
An interesting feature of this system is that for even N 
if £ is an eigenenergy, —£ is an eigenenergy too. Thus, 
the eigcnspectrum exhibits a symmetry, which is different 
from that in the XXX model. The reason is that for even 
N Hamiltonian -ffxY anticommutes with the following 
operator 



Az = 0-12 



'0-3^ 



' C^N-lz- 



(15) 



It is natural to ask if the entanglement and mixedness 
exhibit a similar symmetry. The answer is yes as we will 
see below. 

By exact diagonalization method, we compute the con- 
currence and mixedness numerically in the XY model. In 
Fig. 3, we plot the concurrences for two nearest qubits 
versus eigenenergy per site efe . It is evident that the con- 
currence shows a symmetry with respect to the point of 
ek = 0, which arises due to the fact that the operator Az 
is a local unitary operator, and the entanglement proper- 
ties of the eigenstate with energy £ are the same as those 
of the eigenstate with energy —£. If we restrict ourselves 
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for odd number of qubits, the symmetry breaks. It will 
be interesting to consider entanglement of all eigenstates 
in other physical systems such as quantum chaos sys- 
tem. The study of eigenvalue statistics and eigenvectors 
statistics naturally motivate us to study entanglement 
statistics, which is under consideration. 




FIG. 4: The linear entropy versus eigenenergy per site e^ for 
different A'^ in the XY model. 



to the half part of the eigenspectrum (cfe < 0), the be- 
haviour of the concurrence is similar to that in the XXX 
model, i.e., the concurrence monotonically decreases as 
k increases and there also exist threshold values fcth- 

The linear entropy quantifying the mixedness of the 
two-qubit states is plotted in Fig. 4. It also exhibits a 
symmetry with respect to e/c = 0. And it seems that the 
linear entropy takes its maximum at e^ = for any N 
from Fig. 4 (a-c). However, the partly enlarged version 
(Fig. 4 (d)) of Fig. 4 (c) clearly shows that the maximum 
is not exactly at e^ = for N ~ 10. So, we sec that the 
mixedness reaches its maximum at or near the point of 

In conclusion, we have studied the entanglement and 
mixedness of all eigenstates in two well-established quan- 
tum spin models, i.e., the XXX and XY models. For the 
XXX model, the entanglement and mixedness of a pair 
of nearest-neighbor qubits in a eigenstate are completely 
determined by the corresponding eigenenergy. The en- 
tanglement decreases as we go from ground state to ex- 
cited states, i.e. the more excited the system, the less 
the entanglement. The negativity of eigenenergy implies 
pairwise entanglement, i.e., the entanglement vanishes 
when the eigenenergy is large than zero. In contract 
to the entanglement property, the mixedness exhibits a 
maximum in the middle of energyspectrum, other than 
on the border. From the relation between entanglement 
and eigenenergy, we obtain the finite-size behaviors of 
the entanglement. 

The properties of entanglement and mixedness in the 
XY model are distinct from those in the XXX model. 
For even number of qubits, they both exhibits a sym- 
metry with respect to the zero eigenenergy. However, 



[2; 
[s; 
[4; 
[s; 
[6; 
[7: 

[s; 

[9; 

[lo; 

[11 

[12 
[13 

[14 

[15 

[16 

[17 

[18 

[19 

[20 
[21 
[22' 
[2Z 

[24' 
[25| 



C. H. Bennett, G. Brassard, C. Crcpcau, R. Jozsa, A. 
Peres, and W. K. Wootters, Phys. Rev. Lett. 70, 1895 
(1993). 

T. J. Osborne and M. A. Nielsen, Phys. Rev. A 66, 

032110 (2002). 

M. A. Nielsen, Ph. D thesis. University of Mexico, 1998, 

quant-ph/0011036; 

M. C. Arnesen, S. Bose, and V. Vedral, Phys. Rev. Lett. 

87, 017901 (2001). 

D. Gunlycke, V. M. Kendon, V. Vedral, and S. Bose, 
Phys. Rev. A 64, 042302 (2001). 

X. Wang, Phys. Rev. A 64, 012313 (2001); Phys. Lett. 
A 281, 101 (2001). 

X. Wang, H. Fu, and A. L Solomon, J. Phys. A: Math. 
Gen. 34, 11307(2001); X. Wang and K. M0lmer, Eur. 
Phys. J. D 18, 385(2002). 

B. Q. Jin and V. E. Korepin, Phys. Rev. A 69, 062314 
(2004). 

G. L. Kamta and A. F. Starace, Phys. Rev. Lett. 88, 
107901 (2002). 

K. M. O'Connor and W. K. Wootters, Phys. Rev. A 63, 
0520302 (2001). 

D. A. Meyer and N. R. Wallach, quant-ph/0108104. 
G. K. Brennen and S. S. Bullock, quant-ph/0406064. 
A. Wong and N. Christensen, Phys. Rev. A 63, 044301 
(2001). 

A. Osterloh, L. Amico, G. Falci and R. Fazio, Nature 
416, 608 (2002). 

Y. Sun, Y. G. Chen, and H. Chen, Phys. Rev. A 68, 
044301 (2003). 

J. Vidal, G. Palacios, and R. Mosseri, Phys. Rev. A 69, 
022107 (2004). 

N. Lambert, C. Emary, and T. Brandes, Phys. Rev. Lett. 
92, 073602 (2004). 

G. Vidal, J. L Latorre, E. Rico, and A. Kitaev, Phys. 
Rev. Lett. 90, 227902 (2003) 

S. J. Gu, H. Q. Lin, and Y. Q. Li, Phys. Rev. A 68, 
042330 (2003). 

W. K. Wootters, Phys. Rev. Lett. 80, 2245 (1998). 
X. Wang and P. Zanardi, Phys. Lett. A 301, 1 (2002). 
X. Wang, Phys. Rev. A 66, 044305 (2002). 
G. Jaeger, A. V. Sergienko, B. E. A. Saleh, and M. C. 
Teich, Phys. Rev. A 68, 022318 (2003). 
S. Bose and V. Vedral, Phys. Rev. A 61, 040101 (2000). 
L. Cardy, in Phase transitions and Critical Phenomena, 
edited by C. Domb and J. L. Lebowitz (Academic, Lon- 
don, 1987), Vol. 11. 



